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Abstract: We consider the creation of non-zero Chern-Simons number in a model of the 
early Universe, where the Higgs field experiences a fast quench at the end of inflation and 
subsequently rolls down its potential barrier. Neglecting the expansion, we perform nu¬ 
merical lattice simulations in the Abelian Higgs model in 1-|-1 dimensions with an added 
phenomenological C and P violating term during this stage of so-called tachyonic preheat¬ 
ing. The results suggest that even the sign of the Chern-Simons and thus baryon number 
is dependent on the ratio of the Higgs to W mass. We also discuss the appropriate choice 
of vacuum initial conditions for classical simulations. 


Keywords: Baryogenesis, CP-violation, Out-of-equilibrium field theory. Preheating. 





Contents 


1^ Introduction 

1^ The Abelian-Higgs model 

1^ Classical approximation 

^ Initial conditions 

1^ Numerical results 

|5.1| Dependence on initial conditions 
|5.2| Initial asymmetry 
|5.3| Dependence on mH/mw and n 
|5.4| Suppression of sphaleron wash-out 
|5.5| Volume dependence 
|0| Modelling 

|5.7| Role of topological defects 


s 

I 


M 

H 

0 

p 

0 

0 

P 



1. Introduction 

It was proposed some time ago i § that the observed baryon asymmetry may have 
been produced in the electroweak transition at the end of an inflationary period. The 
mechanism assumed resonant preheating, in which low-momentum bands of the Higgs 
spectrum obtained large occupation numbers. These would then ensure a high effective 
temperature in the low-momentum gauge degrees of freedom, leading to baryon number 
production via Chern-Simons number generation by sphaleron transitions. At a later stage 
the system would thermalize at a low temperature below the electroweak scale, leading 
to the well-known suppression of the same sphaleron processes. Numerical simulations of 
the 1-|-1 dimensional classical Abelian-Higgs model with a C and P violating term in the 
action supported this idea Q. Recently Q, it was proposed that electroweak production 
of Chern-Simons number could even be resonant itself, with a large boosting of the Chern- 
Simons number at preheating. Another scenario that has been put forward exploits the 
production of topological defects in the Higgs field by the Kibble mechanism, resulting in 
a Chern-Simons asymmetry under influence of CP violation Q . 
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We investigate here a related mechanism that does not assume resonant preheating or 
the necessity of topological defect production. Also in this scenario inflation is assumed 
to end at the electroweak scale, leaving the Universe in a cold state. We assume that 
the electroweak transition was sufficiently rapid that it resulted in a spinodal instability 
at essentially zero temperature. In this paper we model it by a quench, in which the 
parameter in the effective Higgs potential +A((/>*(/>)^ changes sign from positive 

to negative on a time scale much shorter than a typical electroweak time. Initially the 
quantum helds are in their semiclassical ground state at > 0. Subsequently /rgg —> — 
and quantum modes in the Higgs field with momenta smaller than /r grow exponentially 
fast, the spinodal instability. After some time large occupation numbers are reached and 
a classical description makes sense. One expects the growing Higgs held to generate a 
growing SU(2) gauge held as well, through the classical equations of motion, and when 
CP-violating interactions are operative, a Chern-Simons asymmetry will be created. Some 
of this may survive or even grow during the ensuing redistribution of energy over the held 
modes, and a practically frozen non-zero Chern-Simons number may emerge if the effective 
temperature is sufficiently low. 

The problem we address ourselves to is: neglecting the expansion of the Universe and 
ignoring remaining inhaton effects, how large a baryon asymmetry is generated in the non¬ 
equilibrium process under inhuence of CP violation? We have in mind the CP violation 
in the Standard Model corresponding to the Cabibbo-Kobayashi-Maskawa matrix, and 
similar terms in an extended model including neutrino mixing. This CP violation is gen¬ 
erally considered much too small, but the arguments leading to this conclusion are based 
on dimensional analysis involving the electroweak symmetry breaking scale or a high 
temperature of order of 100 GeV p, and the situation may be different at zero tempera¬ 
ture and small Higgs condensate during electroweak symmetry breaking. By assuming a 
quenching electroweak transition at zero temperature we may expect to obtain the largest 
possible baryon asymmetry in this kind of scenario. 

In a hrst approach to this problem we study a 1-1-1 dimensional analog model in this 
paper, the abelian-Higgs model with an effective C- and P-violating interaction of the form 
also used in |Q| (for an early study, see |@]) 

Ccv = ( 1 . 1 ) 

The more realistic case of the SU(2)-Higgs model with effective CP-violating term 

£cp = -«</> V TV (1.2) 

will be reported in a future publication (see for preliminary results). Real time numer¬ 
ical studies of the SU(2)-Higgs model in equilibrium including the above CP-violating term 
can be found in Q. An early study of the effect of a quench on the SU(2)-Higgs system 
is in Q. A scenario for the quench is to introduce a scalar held cr, conveniently identihed 
with the inhaton, with a coupling to the Higgs held given by 

V{(j), cr) = cr^ - + A((/>»^ (1.3) 
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i.e. /igg = As the inflaton eventually ends its slow roll and drops to near zero, 

the Higgs symmetry breaking is triggered. This is a hybrid inflation model |10|. It turns 
out that for viability reasons it is necessary to modify it into what is called Inverted Hybrid 
Inflation |Q]; see also for possibilities of inflation ending at the electroweak scale. In 
this study we ignore inflation and the expansion of the Universe. 

The abelian-Higgs model is introduced in section 2. As a consequence of the spinodal 
instability produced by the quench, the Higgs field acquires classical properties, which 
we review in section 3. The basics of this can already be found in [22|. This can be 
exploited to derive realistic initial conditions for numerical simulations using the classical 
approximation, as already advocated in |12]. These initial conditions differ from those e.g. 
in [13|-|2C] in that they do not put power into the short-wavelength modes, which may 
lead to artificial equilibration properties |^, H|. We give here a detailed derivation of 
these initial conditions (section 3 and 4, with a numerical detail in the appendix). Related 
studies are in |^, 241. In section 5 we present our numerical results, give an interpretation 
of these in terms of simple models, and analyze the possible role played by the Kibble 
mechanism. We conclude in section 6. 


2. The Abelian-Higgs model 

The Abelian-Higgs model in 1-|-1 dimensions is given in the continuum by the action 

S = -JdtdxS^ +^- 

+ ( 2 . 1 ) 

where D^cj) = (d^ - iAf,)(l), = (d^ + iA^)(p* , and eoi = +1- We 

have added a P and C breaking term biasing the Chern-Simons number: 

J dxK^e^uF^’''(j)*4> = 27rKNcs4>*4>, ( 2 - 2 ) 

which is the 1-1-1 dimensional analog of the CP-breaking term (|1.2|). The Higgs and W 
masses are given by mjj- = 2iJ?, = (e^/2A) m^. 

We discretize this action on a space-time lattice in the standard way, in a periodic 
spatial volume L. The classical Euler-Lagrange equations of motion for cj), Ai and Aq read 
in temporal gauge (Aq = 0): 


d'odocp = D[Di(j) + (^2 _ 2Xcj)*cj))(P + KdoAicj), (2.3) 

d'doAi = - e^d'oi^cj)), (2.4) 

dJdoAi = e‘^i{(l)*do(l> — do4>*4>) — Kd[{(f>* cj)). (2.5) 

Here is the forward and d'^ is the backward lattice derivative, d^/(x) = [f{x + — 


/(x)]/a^, d(j/(x) = [f{x) — f{x — a^jl)]/a^, with the lattice spacing in the /U-direction. 
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Equation ( p.5[) is the Gauss constraint. In integrated form it enforces zero total charge in 
our periodic volume, 

- docf)*(f)) = 0 . ( 2 . 6 ) 

X 

The Chern-Simons number simplifies to 

Ncs = (2.7) 

X 

It is easy to show that the equations of motion preserve the Gauss constraint and the total 
charge. 


3. Classical approximation 

In the limit of small Higgs and gauge couplings (A/^^ <C 1, <C 1 in 1+1 dimensions), 

we can solve the quantum evolution in the gaussian approximation. Before the quench we 
assume the held to be in its ground state in the potential H((/>) = for simplicity we 

take /To = The quench then just hips the sign of the potential, H(<^) ^ In this 

section we focus on the Higgs held, using the notation appropriate for 1+1 dimensions. 
The generalization to 3+1 dimensions will be obvious. Gauge helds will be included in the 
next section. 

Choosing periodic boundary conditions in space with volume (length) L, we make a 
Fourier decomposition for each real held operator (f>j, j = 1,2, 0 = {(f>i + and 






ikx 


Vl 


(3.1) 


k ^ k 

Modes with wave number will be unstable and grow exponentially, those with 

will not grow. We start by expanding each real held^ in creation and annihilation 
operators just before the quench near t = 0, 

1 

/;;. I ;; I 


:(afce k +a'_^e »= ), tt^ = 


\/2< 




whereas after the quench, for t > 0, we write 

= ake~^‘^k * + (dke^'^k \ = _Aa;“ (dfce"*‘^fc * - *) 

where 


= V±ax2 + P. 


Matching at time zero gives the relation 


1 


Oik = 


/3fc = 


2^ 2ujt 




2^ 2ujt 






1 + ^ dfc + 1 - ^ UL 


OJi. 


Wi, 


1 - ^ dfc + 1 + ^ UL 






W,, 


^We drop the index j. 


(3.2) 


(3.3) 

(3.4) 

(3.5) 

(3.6) 
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For the stable modes is real and /3fc = For the unstable modes is imaginary 

and we write = i\oj^\ (the opposite sign gives equivalent results). Then q;|, = a^k and 
(3\. = /3_fc. In both cases the reality conditions (j)\ = ^-k, = ^-k are satisfied. 

Consider now the unstable modes, i.e. with \k\ strictly smaller than /r and o;^ = 
which grow exponentially fast when S> 1. Neglecting the decaying exponential in 

expressions ( |3.3|) gives 

^k-ake\'^k\\ vTfc « (3.7) 

This strongly suggests classical behavior, since = 0 and consequently (/>fc and 

commute in this approximation. Of course, there are states and observables for which the 
approximation is not valid, e.g. the hermitian operator i[T:k,(j)-i] = 6ki at all times, being 
the canonical commutator. So let us see what happens to the field correlation functions in 
the physically relevant state, the initial state just before t = 0. 

The initial state is assumed to be the ground state just before the quench, |0), which 
satisfies dfc|0) = 0. We can now find the field correlators at time t > 0: 


Cf = {0\Ml\0) 

cr = (OlvTfcvrtlO) 
Ct = Cf* = (0|<^A:Vrl,|0) 




1 + ( ^ - 11 






-2 




u> 


+2 


1 + - 1 COS (w^ t) 


CO, 


COu 


00 


+2 


4io: 


COr 


32-1 sin(2u;. t) + - . 


(3.8) 

(3.9) 
(3.10) 


The correlator (0|7r]!(^fc|0) follows from the commutation relation (^k'^l = '^l^k + i- 

We re-express these correlators in terms of time-dependent particle numbers nk, fre¬ 
quencies ook, and off-diagonal particle numbers hfc, as follows: 


= ink + 1/2)/ook, 

(3.11) 

'k^ — 4" 4/2)iOki 

(3.12) 

= nk + il2. 

(3.13) 


These and oOk have proven to be robust and very useful in numerical studies of interacting 
scalar fields out of equilibrium Often fik is equal to zero, but here it is not and it 

plays an important role in the transition to classical behavior, as we shall see below. We 
note the identity 

[uk + 1/2 nk){nk + 1/2 - hk) = 1/4. (3.14) 

Consider again the unstable modes. Their particle numbers grow exponentially. For 
\oo'^ |t 3> 1 we find 

rik + ^-nk-^^ j;ri^ g2VM2-fc2 4 ^ Wfc « |a;^| = _ f.2^ (315) 
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and the generic field-expectation values behave as classical. To express this more clearly, 
let us introduce sources and which are only nonzero for k in the unstable region, i.e. 
strictly smaller than 

f dx + J^vr) = ^ + J-k^k), Ate = (1 - e)n, (3.16) 


where e is a small positive number <C 1. For sufficiently large times the expectation 
values of products of </>’ and vr’s, with any operator ordering, can be calculated. Using the 
approximation (3.7) we find for the generating functional 


Q[J] = (0| exp 


dx [J^{x)${x) + J’^(x)7r(x)] 


exp 


J dx dy J^{x)Cf{x,y)j'^{y) 


| 0 ) 


(3.17) 


where we summed over a,b = {4>, tt}, and the dominant part of the correlator corresponding 
to the unstable modes is given by^ 




c, 


u k 



k\<lle 


f ^ \^k I ^ _g2\/ fj?-k'^ t 


(3.18) 

(3.19) 


The matrix Cufc is singular, which reflects the fact that irk = \‘^k\4>k in this approxima¬ 
tion. So, for sufficiently large times, the dominant part of the quantum correlators can be 
expressed as a probability distribution in a functional space of classical (f) and vr consisting 
only of unstable modes: 


Q[J] ra M Jldcpdn] 


6{7rk - \<f)k) 


k<fJ.e 


exp 


-\J dxdy(t>^{x)C^lf^{x,y)(l)\y) 


+ J dxJ^{x)ct)^{x) 


(3.20) 


where N is such that ^[0] = 1, and = (j), = tt. The functional measure [dcpdir] can be 

explicitly (and tediously) expressed in terms of the independent real and imaginary parts 
of TTfc and (pk with |A:| < /i^, but we will refrain from doing this here. 


With the help of the classical probability distribution (|3.20D we can calculate expecta¬ 
tion values of products of the usual observables that have a classical correspondence, which 
are represented in the quantum theory by products of (p and vr (symmetrized, if needed, 
e.g. the charge density = —-k^cp’^ + k‘^(p^. Provided such observables are dominated by 
their unstable mode contribution, we can sample the distribution (|3.2C|), use the samples as 


^We used the Campbell-Baker-Haussdorf formula for {0| exp 511 Cfcdfc|0) = exp . 
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initial conditions for subsequent classical dynamical evolution into the non-linear regime, 
and compute expectation values by averaging over the initial conditions. 

The above prescription runs into the question: what to choose for e? It would be nice 
to be able to let e —> 0, but it is clear from ( ^.19 ) that we then have to do something 
different near the boundary of the unstable region, because of the factor \ jsj. We 
propose to make the replacement 


C, 


u k 


ReCk = 


(Uk + ^)/uJk Uk 

flk {Uk + \)uJk 


(3.21) 


which goes over into Cufe for large and to set e = 0. Using the real part corresponds 

to symmetrized products, e.g. ReC^’^ = {Q\\{^k^\ + ^\4>k)\R)i which goes over into the 
corresponding classical correlator in the classical limit |^|. The delta functional in (|3.2C| ) 
is to be omitted in the replacement ( 3.21 ). It is generated automatically in approximate 
form when the instability progresses. This can be seen by introducing the variables 


= J_ 

\/2 VV^ 


^ ( -^= ± (t>k ) , 


(3.22) 


in terms of which the proposed classical distribution takes the form 


exp 


E (—^ 

Z—/ \ -n 1 -\- 1 


+ |2 
k I 


+ 


1^, 


-|2 
k I 


2 \nk + 1/2 + hk nk + 1/2 - hk 


(3.23) 


As time progresses, -|- 1/2 — —> 0 rapidly, for \k\ < /Xe, as illustrated in figure |T|, and 



Figure 1: The difference nk + 1/2 — Uk versus k/fi at early times tfi = 0, 0.3, 0.9, 2.4 (from top to 
bottom at fc = 0). 


the delta functional enforcing oc (vTfc — 0Jk4>k) = 0 in (|3.2C ) appears automatically. This 
squeezing of the (fik — vr^ distribution along Hk = u}k4>k has been studied earlier in detail in 
1^, and in |4|. 
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We could of course drop the modes altogether, since they never grow large. However, 
by including them we will be able to compare with another practical method for obtaining 
initial conditions in classical numerical simulations, which will be introduced below. Here 
we note that the modes with |A:| > /r have + 1/2 ± just oscillating around 1 / 2 , with 
an amplitude decaying like k~‘^. Similarly —> 1/2 for |A;| —> oo. So these modes stay in 
the quantum regime and we do not include them in the initial conditions for the classical 
approximation. Figure ^ shows ± hk and ujk at time t = 4.2 . 



Figure 2: Left: nfc + l/2 + fife (coining down from « 2224 at A: = 0) andnfe + 1/2 —hfc (practically 
zero at fc = 0) versus A:/p, at time = 4.2. Right: jw^Tl and ojk versus A://i at tfi = 4.2. 


We now give an estimate of the time-span for which we have reason to trust the free-field 
approximation. Consider the unstable mode contribution to (0|i/)^(x)|0), (0|(/^(x)|0)unst = 
if?'. We define a time tni (‘t-^non-linear’) such that y? has grown so large that it is at the 
inflection point of the potential, i.e. d‘^V[ip)j = 0, or p?' = /i^/(3A) = u^/3. For times 
t > tni) non-linearities will certainly come into play. Using the correlator ( |3.8| ) in the 
infinite volume limit this gives the criterion, in d spatial dimensions. 




\k\<^,{27:Y ^ 3A- 


(3.24) 


Figure |3| (left) shows as a function of time. For large time it behaves as 




Figure 3: Left: the function p'^{t )/versus tp for d =1 (full), 3 (dashed). It reaches the value 
3 at t/r « 2.05, 4.20, respectively for d = 1, 3. Right: plot of n/^^(2.05) (full) and of n^'''^(4.2)/100 
(dashed) versus k/p. 

Cd , with Cl = 0.070528, C 3 = 0.0056111, which is accurate within 5% 
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for /if > 2, 2.5 in one, respectively three dimensions. For the Standard Model, A = 
Ri 1/9 for m\/m^ ^ 2, and the corresponding value 3 « 1/(3A) is indi¬ 
cated in the figure by the horizontal line. We see that in 1-|-1 dimensions this line is reached 
earlier. Figure ^ (right) shows the particle numbers at the time fnij for d = 1 and 3. 
For small k they are much larger than one, especially in three dimensions, which justifies 
switching to the classical approximation at a time somewhat before fni- 


4. Initial conditions 


Our basic assumption is that at the end of inflation and before the instability has set 
in, the system is in its time-dependent semi-classical ground state, which we approximate 
by the free vacuum of the Higgs and gauge fields corresponding to /r^g > 0. As the 
unstable Higgs fields grow large, we switch to the classical description at a roll-off time fro 
before non-linearities have become important, fro < fni- This means that we sample the 
distribution ( 3.23|) , construct (j){x) and 7r(x) from the and use these as initial conditions 
for subsequent classical evolution. Since all Higgs field modes that have not grown large 
are neglected this way (those with \k\ > /x), it is natural to set the initial gauge field to 
zero as well (i.e. Ai^ = Aik = 0 for all k). However, this would wrongly ignore the Gauss 
constraint. 

The global Gauss constraint (|2.6| ) of zero total charge renders the distribution ( 3.23| ) 
non-gaussian. We take this into account by Monte Garlo methods. See the appendix for 
more details of this. Given a sample of the initial (j) and vr with zero total charge, we 
then keep Ai = 0 and satisfy Gauss’ law by solving for the non-zero modes of Ai (its zero 
momentum mode is set to zero as well). Note that this brings the coupling into the 
initial conditions. 

For f > fro the classical evolution takes place with non-zero and A. But if the 
non-linear interactions are small at the time fro, then at first the correlators and 

are still given by the expressions (^(^ ) ~ (|3.10|) . So the results should not depend on the 
precise value of fro- In fact, since the classical evolution for (7“^ is identical to the quantum 
evolution in the gaussian case, we may as well send fro to zero. This might give better 
results in case of stronger couplings, for which the non-linearities cannot be neglected even 
at small times. Since = hfc = 0 (and ojk = w/") at f = 0, the denominators in the 
distribution ( [1.23 ) only contain the factor 1/2 — which is why we call this choice of initial 
conditions the ‘just a half’ method. This is almost identical to the choice made in [p!^- |p0|| , 
but with the important difference that we have a cut-off on the initialized modes that are 
non-zero: \k\ < fi. 

A completely different initialization is obtained by replacing -|- 1/2 and in ( 3.23| ) 
by a Bose-Einstein (BE) distribution. 


+ 1/2 


1 




hk 0 , 


(4.1) 


without restriction on k. Such an ensemble of quantal BE initial conditions for a classical 
approximation may seem strange, but they do indeed represent a free-held thermal quantum 
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density operator |^|, i.e. of the form 


p oc exp 



(4.2) 


for each real component of the Higgs field. The temperature T controls the initial fluctua¬ 
tions that start up the tachyonic instability. Choosing T low enough, the high-momentum 
modes are sufficiently suppressed to avoid Rayleigh-Jeans problems or regularization arte¬ 
facts, at least for some time after start-up, since classical equilibration to equipartition is 
a slow process. 

In summary, we have considered three types of initial conditions for the classical evo¬ 
lution: 


‘spinodal’: sampling (p.23) at roll-off time 2/r ^ < fro < ^ni (assuming > 2p ^), 


‘just the half’: sampling ( 3.23 ) with Uk = hk = 0 at t = 0, 
thermal: sampling (3.231) with the BE form ([4.1j), with T//r <C 1. 


5. Numerical results 


Most of the results presented here are for a number of lattice points N = 512, lattice 
spacing amn = ay/2p. = 0.3, where a = ai, temporal lattice spacing oq = 0.1 oi, and 
volume Lmu = 153.6. We have typically used an ensemble of 1000 initial conditions and 
have run for a time t = 600 keeping track of the average (</>^</>) and Chern-Simons 
number, (Ncs)- We studied the dependence on the coupling ratio 2A/e^ = 
and various coupling strengths = 1/8, ..., 1/512. We also checked, using larger 

lattices {N = 10240) that the final asymmetry is proportional to the physical volume, 
{Ncs) oc mnL. 


5.1 Dependence on initial conditions 


Figure 1^ shows and (iVcs) at early times, for spinodal initial conditions with various 

choices of the roll-off time tj-o, and a comparison is made with ‘just a half’ and thermal 
initial conditions. The CP-asymmetry parameter k = —0.03. The coupling is fairly weak, 
p^/X = 8, for which fni « 1.98 (from ( 3.24 ) for ifl = 8/3), with 2A/e^ = 1 {mu = mw)- 
We see the curves for ((/>^(/) already breaking away from the ‘just a half’ curve for somewhat 
smaller values of fro than tni- We have checked that these deviations diminish for weaker 
couplings. The plot for (iVcs) contains also a comparison with the thermal method for 
T/niH = 0.1 and 0.3. Compared to the ‘just a half’ curve (tro = 0), the thermal curve has 
a longer ‘waiting time’ before the initial dip occurs, which indicates smaller fluctuations in 
the initial conditions. Apparently, the effect of the vacuum fluctuations (characterized by 
Uk + 1/2 = 1/2) is substantial in comparison with at T/mu = 0.1. 
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Figure 4: Top: Dependence of on the roll-off time fro with ‘spinodal’ initial conditions, 
including ‘just a half’ {t^o = 0). The ‘thermal’ case with T/mn = 0.1 is also shown. Bottom: {Ncs) 
in runs with three different initial condition schemes. 

5.2 Initial asymmetry 

We can estimate the initial effect of the C and P breaking term from the equations of 
motion. Approximating the Higgs field by its homogeneous mode falling in the inverted 
quadratic potential, = (/>(0) exp(^t), inserting this into the equation of motion 
for Ai, with Ai(0) = Ai(0) = 0, neglecting the current-term (—2e^|(/>pAi) and taking into 
account only the «:-term, we find at t 

Ai{tni = -e‘^ [ (l>*c/)dt ^ (5.1) 

Jo 

where we have neglected the seed |</>(0)p compared to |</>(tni)P- Using cf) = ((/>! -|- i(f) 2 )/V^ 
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and choosing the 1-axis in the complex plane along (j) (i.e. (j )2 = 0), gives |</>(tni)P = V^i/2 = 
IJ?/6X (cf. ( |3.24 )), and the estimate 

(5,2) 


iVcs{t„,) = f — 
2tt \mw 


6\/2 ■ 


The actual value of the asymmetry in the full simulation in the minimum of the initial 
dip in Nqs in figure is indeed within a factor of two of this estimate (figure , and in 
particular, the sign is right. 



Figure 5: Examples of (iVcs) for mn/mw = 0.625 (top) and 1.0625 (bottom). The volume 
LrriH = 153.6 and k = —0.05. 


5.3 Dependence on mH/rriw and k 

After the initial dip, the behaviour of the Chern-Simons number depends on the ratio of 
the Higgs to W mass, as shown in hgure |^. For niH/niw around 0.625 (Nqs) keeps the 


- 12 


















same sign as the initial dip (top plot), for most others it ends up with the opposite sign 
(bottom plot). The final Nqs is typically larger than the value at the bottom of the initial 
dip due to what looks like resonant behavior. 

We end the simulation at a time mni = 600 when the trajectories are stuck (see below), 
and measure the final value of the average Chern-Simons number, (Ws)- The distribution 
of the 1000 trajectories in such an average is just gaussian (figure P). 



Ncs 

Figure 6: The distribution of the Chern-Simons numbers at irnn = 600 in an ensemble of 1000 
initial conditions, k = —0.05, mH/mw = 1-0. 

The dependence of the final Chern-Simons number on the mass ratio niH/mw at fixed 
ttihL and k (so changing e^) is quite complicated, see figure 0. Only in a narrow 
range around niH/mw = 0.6 is the sign of (Nqs) the same as that of the initial dip and 
the input k (negative). For large mass ratio the asymmetry is expected to go down because 
of the explicit factor accompanying k in (^.4|). However, the strong dependence on the 
mass ratio is related to the resonant behavior mentioned previously. 

The actual value of k in realistic models is presumably quite small, so it is comforting 
to see that the behavior of the final Nqs as a function of k is simply linear, see figure ^ 

5.4 Suppression of sphaleron wash-out 

The effective temperature at the end of tachyonic preheating should be low enough that 
any Chern-Simons asymmetry created by the tachyonic instablility does not get washed 
out through equilibrium-type sphaleron processes. If the system ends up deep enough in 
the broken phase the sphaleron rate will be exponentially suppressed with temperature. 
We can estimate the rate by running single configurations for a long time and monitor 
effective temperature and transition rate. Figure ^ shows some examples. 

The Chern-Simons number eventually gets almost stuck at times of order 1000 
and afterwards the transition rate is strongly suppressed. We estimate the temperature of 
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Figure 7: Dependence of the final (iVcs) on mH/mw- 



Figure 8: Dependence of the final (Ncs) on k for mn/mw = 1- 


the system by monitoring the average canonical momentum of the Chern-Simons degree 
of freedom (this is the only dynamical degree of freedom in the gauge field, as 

is clear e.g. from the formulation in the Coulomb gauge [p^]). A typical evolution of this 
temperature is shown in Figure 10, For a somewhat different value of niH/mw h was shown 
in 129] that the sphaleron rate in this model is exponentially suppressed by the sphaleron 
Boltzmann factor when (3' = /{XT) is larger than ~ 7. In the case of figure |I^ we get 

(3' ~ 10 and so we are deep in the broken phase with a strongly suppressed sphaleron rate. 
This is also clear from estimates of the rate from pictures like figure |9[ When averaging over 
initial conditions, it turns out that the average Chern-Simons number gets stuck already at 


- 14 











Figure 9: Examples of Ncs in single trajectories, mn/raw = 1- 



Figure 10: The effective temperature calculated from the canonical momentum of the Chern- 
Simons degree of freedom. 

tmu ~ 200 (fig. 0)- Its distribution (fig. then widens at a rate equal to the sphaleron 
rate. 

Of course, in a realistic application the additional degrees of freedom in the Standard 
Model and the expansion of the universe are expected to lead to further suppression of the 
rate, such that it ends up being negligible at practically zero temperature. 

5.5 Volume dependence 

When the volume is increased, the fluctuations of the final Chern-Simons number over 
the initial ensemble also grow, but the fluctuations in the density Nqs/L should go down 
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Figure 11: (Ncs) for different k at mn/mw = 1, Lmn = 153.6, except for the top curve, which 
represents iVcs/20 in the 20 times larger volume Lrufj = 3072. 


like . For a very large volume, one classical realization should suffice to provide an 

accurate estimate of the density (which in the realistic case would determine the fermion- 
number density of the universe). We have performed simulations using a 20 times larger 
volume {N = 10240, Lmn = 3072), which should result in an average density that is 
approximately the same as presented in the previous sections (using Lmn = 153.6) , with 
a standard deviation smaller by a factor y/W. We have checked this by using also 20 times 
fewer initial conditions (1000 ^ 50), and found the same standard deviation for the density 
as before. See the upper two curves and their final-time error bars in figure 11 (the top 
curve represents the largest volume). A closer look revealed that (up to an insignificant 
shift in time) the two curves are actually indistinguishable until tniH ~ 13. Also, the 
lattice spacing dependence is very small. 




Figure 12: Analog of figure || in the homogeneous approximation; dashed: (jp' versus imn, con¬ 
tinuous: Ncs- Left: niH/mw = 0.625, right: mH/mw = 1.0625. The initial values are </> = 0.02, 
(j) = Ai = Ai = 0. 
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5.6 Modelling 


In this section we attempt to interpret the data in terms of simple models for the dynamics 
of Ncs, or equivalently, the homogeneous gauge field in the Coulomb gauge, Ai. The first 
model is simply the restriction of the equations of motion to homogeneous fields, with the 
full non-linear dynamics. Without loss of generality, we can assume (p to be real. Figure 


the non-linearity of the oscillations in the fields, due to the absence of damping by the in¬ 
homogeneous modes. However, it appears that the behavior until the second minumum of 
{4>*4>) {truH ~ 12) is reasonably well represented. The initial dips in figure |5| ([T^) are -0.66 
(-0.66) and -0.74 (-0.53), respectively for niH/mw = 0.625 and 1.0625. In the former case 
the Chern-Simons number oscillates more rapidly when (?!*(/) is large, because the W mass 
is larger, and when (p*^ is low again after its first maximum, A^cs coasts along almost freely 
in the negative (positive) direction in the left (right) plot of figure |T^. The homogeneous 
approximation breaks clearly down already in the region of the second minimum of 
but the sign of the true final asymmetry (fig. ^) is the same as the sign of iVcs in this 
region (12 < tnin ^ 20) in the homogeneous approximation (fig. |^. 

When we add damping terms jp and 7 a^i to the homogeneous approximation, the 
resulting can be made to look pretty much like figure including the oscillation period 
. However, the resulting A^cs then simply performes a damped oscillation around zero. In 
the real simulations gets stuck in a minimum between the sphaleron energy barriers. 
We have tried to model this by adding a periodic potential Vs{LAi,p), periodic in LAi 
with period 27r, and height equal to the sphaleron energy Eg = (2m///3)u^, where is the 
expectation value 2p*p = in the classical ground state. The effective action 

then takes the form 


12 shows the result for parameter values as in figure g. A striking difference with fig. ^ is 


•S'eff = L dt 


2e2 


p fJ- p — \p -\-i^A\p -\-V{LAi, p') 


(5.3) 


Since the homogeneous p is varying in time we try the simplistic form V{LAi, p) oc 


V{LAuP) = P^L-^f{LAi), 


(5.4) 


with the tentative conditions 

f{LAi) = ^mnL, LAi = vr, (5.5) 

= {LAif+ 0{{LAif), LAi^O. (5.6) 

The first condition represents the sphaleron energy, the second takes care of the fact that 
the effective Ai mass is 2e^p‘^. A possible solution can be given in the form f{LAi) = 
Yl'li=o^n{'fnHL) cos(nLAi), with coefficients Cn depending on mnL. The equations of 
motion follow straightforwardly from the effective action, after adding also the required 
damping terms, 

p + 'yp = (mHl2 - 2\p^ -g 2L“^/(27rA'cs) + P, (5.7) 

Ncs + le^Ncs = -(e^/2vr) (^f'{2 'kNcs)P‘^ + 2kLPP^ (5.8) 
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where we assumed the damping for Ai to be proportional to e^. 

However, for muL of order 1 (the sphaleron size) and n = 0.05, the K-term is much 
too small to push the Chern-Simons number over the sphaleron barrier and the resulting 
asymmetry is zero. For larger volumes the barrier itself is too high to get A^cs sufficiently 
away from zero; the /'-term in ( ^. 8 | ) becomes very large near the top of the barrier (cf. 
(^)). The /-term in ( |5.7|) may actually be neglected for large L. 

To proceed, we shall interpret the above equations as effective equations for the Chern- 
Simons density, with an effective potential / that is not constrained by the barrier condition 


(|5.5|) , and that is even time-dependent. We have to incorporate the fact that at early times 
t < the quadratic form L~‘^f{LAi) = A\ gives a reasonable description of the 

data. The bottom of the dips in figure |T^ is deeper than the sphaleron value —1/2. This 
means that the sphalerons have not appeared yet, because A^cs bounces back (on the 
potential) instead of rolling down the hill on the other side of a (lower) sphaleron barrier. 
We model this by an effective potential that changes from quadratic to periodic, e.g. 


L 


-2 . 


f(LAi) r^[2 - 2cos(Mi)], (5.9) 

i{t) = + ^\[l + tscnh{5{t — tQ))]/2, (5.10) 

with parameters such that l{t) is very small at t = 0 (and consequently L~^f{LAi) ~ 
and reasonably large at f = oo, such that Ai is able to hop over the barrier. For simplicity 
we set = 153.6 the value of L in most figures. Then, with = 0.01 to = 
5 = 0.16 m//, and damping coefficients 7 = 0.14 mj/, 7 = 0.8 we get the 
result shown in figure 1^. The other parameters are as in figure ^ and 12. 
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Figure 13: Analog of figure including the time-dependent sphaleron-like potential 
dashed: versus tru/r, continuous: JVcs- Left: niH/mw = 0.625, right: niH/rnw = 1.0625. 

Although figure |l^ appears to capture the qualitative behavior of the process, we should 
keep in mind that the model system is quite chaotic, and a small change of e.g. damping 
coefficients can change the result. This could be avoided by considering an ensemble of 
initial conditions. The real simulation lacks of course the arbitrary parameters of the 
modelling, and it is also more subtle, e.g. in that the resonance-like behavior seen in figure 
1 ^ is not very well captured by the modelling. 

5.7 Role of topological defects 

The mechanism proposed in |^, ^ assumed that the electroweak transition produced topo- 


(5.9,5.10); 
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logical defects in the Higgs field, by the Kibble mechanism, with winding numbers averaging 
to zero in absence of CP violation. With CP-violating there would then subsequently be a 
bias towards non-zero average winding number, with a corresponding Chern-Simons num¬ 
ber characterizing the vacuum. Translated to our current notation the asymmetry would 
then be given by a formula of the form ncs = ^cs/L = CKUd, where Ud is the density of 
defects and c a dimensionless factor. 

However, although we do not doubt the fact that there will be a non-zero density of 
topological defects shortly after the transition, this does not play a role in our interpretation 
of the asymmetry as presented in the previous section. With the effective C and P violation 
used here, the all important initial asymmetry (section |5.2D is already about 10% or more 
of the final value, and it has nothing to do with the winding-number density of the Higgs 
field; it depends only on kcIq (</>*</>) ■ 

200 

150 

100 

50 

0 

0 10 20 30 40 50 60 0 10 20 30 40 50 60 

nij^t 

Figure 14: Winding number of the Higgs field Nw, Chern-Simons number Ncs and in two 
large-volume [LmH = 3072) simulations, one with the usual thermal initial conditions (left), and 
one for which the initial Higgs winding is suppressed (right). 

As a check we performed a simulation with initial conditions such, that there is no 
Kibble mechanism. They were obtained from the thermal initial values by rotating (p{x) 
and 7r(x) to the positive real axis, (i.e. (Re</>,Imiji)) ^ (|iji|,0), (Re7r,Im7r)) —> (|7r|,0)), 
and subsequently enforcing Gauss’s law. This has the effect that initially there is no 
topological winding number in the Higgs field. Figure |T^ shows an example of the effect on 
one large-volume configuration, for the mu = mnw case. We have also plotted the winding 
number of the Higgs field^. 

+ ai) - a{x)Y (5.11) 

ZTT 

X 

with (/>(x) = p{x)e'^°‘^^\ and where [.. denotes taking the value modulo 27r. The artifi¬ 
cially non-winding initial configurations are evidently much less random, and we see indeed 
that (/>*(/) in the right-hand plot of fig. 0 suffers less damping in its first half-period than 
in the left-hand plot (which resembles closely the corresponding plot in fig. ||). Note also 

®Using the gauge-invariant winding number of gave identical results 
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the deeper and broader minimum around imn = 13, which reminds us of that in figure 12 
However, we note in particular the complete absence of Higgs-winding in the right-hand 


plot of fig. 14 in 0 < tmu < 12, as compared to the somehat noisy winding in this time- 
interval in the left-hand plot. After time tniH — 30, it appears that the Higgs-winding 
number is pulled along by the Chern-Simons number (bringing down the magnitude of 
the covariant derivative Di4>), such that in both plots Nqs ~ Despite the large ini¬ 
tial deviations, the final asymmetry is semi-quantitatively unchanged with the artificially 
winding-suppressing initial conditions. 

We conclude that it is the gauge field that, after being biased into the ‘initial dip’, pulls 
the Higgs phase along, and the picture of a multitude of topological defects un-winding 
under a biasing C(P) asymmetry in the equations of motion is not relevant for the final 
asymmetry in the present model. 


6. Conclusion 


The H-l D abelian-Higgs model illustrates nicely that a sizable Chern-Simons asymmetry 
can be produced by a tachyonic electroweak transition under the influence of the usual 
effective C- and P-violating interaction. The numerical results can be summarized as 


(Ncs) r. ^7 
ncs = ^ = -0.7 ktuh, 


( 6 . 1 ) 


at mn /rnw = 1 with the dependence on Higgs to W mass ratio as shown in figure The 
linearity in k enabled us to carry out the simulations at a much larger value than might 
be expected in realistic applications. 

Besides initial conditions that are motivated by the quantum-to-classical transition in 
the gaussian approximation, we also used low-temperature thermal noise for generating 
initial configurations. We found that the quantitative results do depend moderately on the 
choice of initial conditions, but not the qualitative outcome. 

It came as a surprise to us that (even the sign of) the final Chern-Simons number and 
thus the corresponding baryon number asymmetry is very sensitive to the Higgs mass. An 


interpretation of this intriguing effect was given in section 5.6 


The mechanism for the generation of the asymmetry here is different from that sug¬ 
gested in ll], 1^ 1^ , since neither resonant preheating nor Kibble-like generation of topological 
defects plays a crucial role in the model studied here. 

We have neglected here any effects related to the dynamics of the expansion of the 
universe, because our primary aim is to see the order of magnitude of the asymmetry that 
can be generated, given a form of CP violation. It will be very interesting to see similar 


results in the physically relevant SU(2)-Higgs theory in 3-1-1 dimension |25|. 


A. Implementing zero total charge 


As discussed in section we need to generate the distribution ( 3.23|) subject to the con¬ 
straint of zero total charge, e^((/>*7r — 7r*4>), or in terms of Fourier variables (defined as 
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in (|3l|) ), 




- vr: 


2 * i,l 


= 0 . 


(A.l) 


On the spatial lattice with an even number of N sites, x = ma, m = 0, ..., N — 1, the 
wave vectors can be chosen to take the values 


k = 


27rn 


n = -N/2 + l,...,N/2. 


(A.2) 


The reality of the fields and tt^{x), j = 1,2, and the fact that ex.p{ikx) is real for 

n = 0, N/2, imply that we can write 



(A.3) 

vrifc = (4 + * 4 ) > n = 1,..., y - 1, 

(A.4) 

yy 4 , n = 4 , n = 0, 

(A.5) 




where the real a’s, ..., c’s, are independent variables. In terms of these the zero-charge 
condition (A.l) takes the form 


Ar/2-l 


Q — Co«0 “ CqOo -I- clfj2a1fj2 — Cw/2®]v/2 + + ^k^k ~ “ 0. (A.6) 


k=l 


To implement the zero charge constraint, the probability distribution (3.22) has to be 



Figure 15: The distribution of (plain) (/)J (the constrained variable) and of (checkered) ^ 

of 500 initial conditions for a thermal ensemble at temperature T/mn = 0.1. Overlaid, a gaussian 
with the relevant width. 
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multiplied by 6{Q{a,b,c,d)), which means that it no longer depends quadratically on the 
aj^,. ■ ■, cj^- For example, for thermal initial conditions hk = 0, and we have 


P{a,b,c,d) (X 6{Q)ex.p 


-E 

2 ^ 


no + 1/2 njy/2 + 1/2 j 


(A.7) 


We can now integrate out one variable, say aj to get rid of the 5 function 




2 ^ 

k,j 

^ nk + 1/2 J 


1 V 

2 yriQ + 1/2 n-N/2 + 1/2 y 2 ng + 1/2 


(A.8) 


where 

Q(a, b, c, d) = Q{a, b, c, d) + Coag. (A.9) 

This distribution is no longer gaussian in the remaining a^, ..., c^. We sample it using 
Monte-Carlo methods. Sampling a distribution with a 5-function is notoriously difficult. 
However if we have enough variables on which to “distribute” the constraint, the deviation 
from a gaussian distribution is expected to be small, which is indeed the case, see figure 

an ideal Monte-Carlo algorithm. As it turned out, the number of variables was sufficiently 
large that no problem was encountered. 

Having produced a realization, we solve for Oq and so construct a field configuration 
{4d (x), TT^ (x)} with zero Q. As mentioned in section]^, the gauge field configuration then 
follows from Ai(a:) = 0 and determining (9gAi(a:) by solving the Gauss constraint. 


15. Our choice of dependent variable (oq) should not matter, in principle, but only with 
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